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Abstract 

Over a field or skew field F with an involution a (possibly the 
identity involution), each singular square matrix A is *congruent to a 
direct sum 

5**^5 = J„, ©■•■©J„p, l<ni<---<np, 

in which S is nonsingular and S* = S'^; B is nonsingular and is deter- 
mined by A up to *congruence; and the n^-by-ni singular Jordan blocks 
and their multiplicities are uniquely determined by A. We give a 
regularization algorithm that needs only elementary row operations 
to construct such a decomposition. If F = C (respectively, F = M), 
we exhibit a regularization algorithm that uses only unitary (respec- 
tively, real orthogonal) transformations and a reduced form that can 
be achieved via a unitary ^congruence or congruence (respectively, a 
real orthogonal congruence). The selfadjoint matrix pencil A + \A* is 
decomposed by our regularization algorithm into the direct sum 

S*{A + \A*)S = {B + XB*) © {J„, -h AJ*J © • • ■ © (J„^ -I- AJ*J 

with selfdajoint summands. 
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1 Introduction 



All of the matrices that we consider are over a field or skew field F with an 
involution a i— > a, that is, a bijection on F such that 

a + b = a + b, ab = ba, a = a. 

If F is a field, the identity mapping a i-^ a on F is always an involution; over 
the complex field, complex conjugation a i— > a is an involution. We refer to 
a as the conjugate of a. 

The entry- wise conjugate of the transpose of a matrix A = [aij] is de- 
noted by 

If there is a square nonsingular matrix S such that S*AS = B, then A and B 
are said to be *congruent; if the involution on F is the identity, i.e., S* = S'^ 
and S*AS = S^AS = B, we say that A and B are congruent. Congruence 
of matrices (sometimes called T-congruence) is therefore a special type of 
*congruence in which the involution is the identity. Over the complex field 
with complex conjugation as the involution, *congruence is sometimes called 
conjunctivity. If A is nonsingular, we write A~* = {A*)~^. 
Let 

1 





denote the n-by-n singular Jordan block. 

For any m-by-n matrix A (that is, A € F™^") we write N{A) := {x G 
F" : Ax = 0} (the null space of A) and denote its dimension by diniN(A) = 
nullity j4. If A is square, we let 

^[''1 := A • • • e A {k times). 

In Section [2] we describe a constructive regularization algorithm that 
determines a regularizing decomposition 

B © Jni ® • • • © Jup, B nonsingular and 1 < ni < ■ ■ ■ < Up (1) 

to which a given square singular matrix A is *congruent. The *congruence 
class of B (the regular part of A under *congruence) and the sizes and mul- 
tiplicities of the direct summands Jm , • • • , Jup (the singular part of A under 
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*congruence) are all uniquely determined by the *congruence class of A. If 
F = C (respectively, F = M), the regularizing decomposition ([1]) can be de- 
termined using only unitary (respectively, real orthogonal) transformations. 
Our proof of the existence and uniqueness of the regularizing decomposition 
([!]) uses two geometric *congruence invariants that we discuss in Section O 
dimiV(A) and dim{N{A*) n N{A)). 

In Section m we exhibit a canonical sparse form that is *congruent to A 
and determines the sizes and multiplicities of the nilpotent direct summands 
in the regularizing decomposition ([T]). The essential parameters of the sparse 
form are identical to those produced by our regularization algorithm, which 
verifies the validity of the algorithm. When F = C or M, we describe a 
reduced form related to the canonical sparse form that can be achieved 
using only unitary *congruences or T-congruences. 

The regularization algorithm reduces the problem of determining a *con- 
gruence canonical form to the nonsingular case. A complete set of *congru- 
ence canonical forms (up to classification of Hermitian forms) when F is 
a field with characteristic not equal to two is given in [Sj Theorem 3]; see 
also [U Theorem 2]. A nonalgorithmic reduction to the nonsingular case 
was given by Gabriel for bilinear forms [3]; his method was extended in f7\ 
to sesquilinear forms, and in to systems of sesquilinear forms and linear 
mappings. The form of the regularizing decomposition ([T|) is implicit in the 
statement of Proposition 3.1 in [ij when F is a field and the involution is the 
identity; the construction employed in its proof does not suggest a simple 
algorithm for identifying the parameters in ([T|). 

If A,B G fnT-xn^ then the polynomial matrix A + XB is called a matrix 
pencil. Two matrix pencils A + XB and A' + XB' are said to be strictly 
equivalent if there exist nonsingular matrices S and R such that S{A + 
XB)R = A' + XB'. Van Dooren ^lOj has given an algorithm that uses only 
unitary transformations and for each complex matrix pencil A+XB produces 
a strictly equivalent pencil 

{C + XD) e (Ml + XNi) e • • • © (M; + XNi) (2) 

in which C and D are nonsingular constituents of the regular part C + XD 
of the Kronecker canonical form of A + XB; each Mj + XNi is a singular 
direct summand of that canonical form [U Section XII, Theorem 5]. Each 
Mi + XNi has the form 

In + XJn, Jn + XIn, Fn + XGn, Or + XF^ 
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for some n, in which 



1 







1 



and Gn 



1 







1 



are (n — l)-by-n. 



The direct sum ([2]) is a regularizing decomposition oi A + \B; C + XD is the 
regular part oi A + XB. Van Dooren's algorithm was extended to cycles of 
linear mappings with arbitrary orientation of arrows in 

If Van Dooren's algorithm is used to construct a regularizing decomposi- 
tion of a *selfadjoint matrix pencil A + XA*, the regular part produced need 
not be *selfadjoint. However, the regularizing decomposition of vl + XA* 
that we describe in Section [5] always produces a *selfadjoint regular part. 

For any nonnegative integers m and n, we denote the m-by-n zero matrix 
by Omn, or by Om if = n-. The n- by-zero matrix 0„o is understood to 
represent the linear mapping ^ F"; the zero-by-n matrix Oon represents 
the linear mapping — >■ 0; the zero-by-zero matrix Oq represents the linear 
mapping — > 0. For every p x q matrix Mpg we have 



pq 







mO 



'Mpg 


■ 




'Mpg 


OpO 




'Mpg 





OmO. 






OmO. 







and 

Mpg 

In particular, 



' Oon 



'Mpg 


" 




'Mpg 







Oon 




_Oog 


Oon 







'pO 



Ooq = 



'pq 



pq OpnJ 



and Jn = Oq • Consistent with the definition of singularity, our convention 
is that a zero-by-zero matrix is nonsingular. 



2 The regularization algorithm 

The first stage in our regularization algorithm for a singular square matrix 
A is to reduce it by *congruence transformations in two steps that construct 
a smaller matrix j4(i) and integers mi and m2 as follows: 

Step 1 Choose a nonsingular S such that the top rows of SA are linearly 
independent and the bottom mi rows are zero, then form {SA)S* and 
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partition it so that the upper left block is square: 



'A'' 




SA = 
I — > SAS* 
The integer mi is the nullity of A. 



{S is nonsingular and the rows of 
A' are linearly independent) 

{S is the same and 



'A'S* 




" M 


N 













M is square) 



(3) 



Step 2 Choose a nonsingular R such that the top rows of RN are zero and 
the bottom m2 rows are linearly independent: 



RN 




E 



{R is nonsingular and the rows of 
E are linearly independent) 



(4) 



The integer m2 is the rank of N . Now perform a *congruence of S*AS 

(5) 
(6) 



with R © 


/: 




" M 













{R®I) 


" M 


iV 








RMR* 


RN 













(i?©/)* 






"Ad ^ 







C7 D 




}m2 







}mi 



The block RMR* has been partitioned so that D is m2-by-m2. The 
size of the square matrix A(^i'j is strictly less than that of A. 

If ^(1) is nonsingular, the algorithm terminates. If ^4(1) is singular, the 
second stage of the regularization algorithm is to perform the two *congru- 
ences ^ and ([5]) on it and obtain integers (the nullity of ^(i)) and 772.4, 
and a square matrix A^2) whose size is strictly less than that of ^(i)- 

The regularization algorithm proceeds from stage k to stage A; + 1 by 
performing the two *congruences ([3j) and ([5]) on the singular square matrix 
^(fc-i) to obtain m2k-i, rn2k-, and ^(fc). When the algorithm terminates at 
stage r with a square matrix Ai^^^ that is nonsingular, we have in hand a 
non- increasing sequence of integers mi > 7712 > • • • > m2T-i > m2T > and 
a nonsingular matrix ^(r)- Our main result is that these data determine the 
singular part of A under *congruence as well as the *congruence class of the 
regular part according to the following rule: 
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Theorem 1. Let A he a given square singular matrix over F and apply the 
regularization algorithm to it. Then A is *congruent to A(^^^ M , in which 
j4(t-) is nonsingular and 

JYf — j[™i-™2] ^ j[m2-m3] ^ jlm3-m4] . . . ^ j[m2r-i-m2T] ^ j[m2T] /^n 

The integers rui > m2 > • • • > ?712t-i > n^2T ^ 0, as well as the *congruence 
class o/yl(T-), are uniquely determined by the *congruence class of A. 

In the next section we offer a geometric interpretation for the integers 
rui in ([7]) and explain why they and the *congruence class of each of the 
square matrices ^(^^ produced by the regularization algorithm are *congru- 
ence invariants of A. Implicit in the regularization algorithm are certain 
reductions of A by *congruences that we refine in order to explain why the 
regularizing decomposition in ^ is valid. 

The nonsingular matrices S and R in the two *congruence steps of the 
regularization algorithm can always be constructed with elementary row 
operations. For the complex (respectively, real) field, it can be useful for 
numerical implementation to know that S and R may be chosen to be unitary 
(respectively, real orthogonal). 

Theorem 2. Let A be a given square singular complex (respectively, real) 
matrix. The regularizing decomposition ^ of A can be determined using 
only unitary (respectively, real orthogonal) transformations. 

Proof, (a) Suppose F = C with complex conjugation as the involution. Let 
A = U*TjZ be a singular value decomposition in which S = Si 0^,^, Si 
is positive diagonal, and U and Z are unitary. The choice S = U achieves 
the required reduction in Step 1. In Step 2, let N = V*T,W be a singular 
value decomposition in which V and W are unitary, S = Si 0, and Si is 
positive diagonal and m2-by-m2. Let 



1 



1 



be the reversal matrix whose size is the same as that of V. Then N = 
{PVy{Pt)W, (PS)VF has the block form dH), and V := PV is unitary, so 
we may take R = V m. Step 2. Thus, A is unitarily *congruent (unitarily 
similar) to a block matrix of the form ([6]) in which D is square and each of 
E and [^(i) B] has linearly independent rows. 
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(b) Suppose F = C with the identity involution. In Step 1, choose S = U 
from (a). In Step 2, choose R = V from (a). Thus, A is unitarily T- 
congruent to a block matrix of the form ([6]) in which D is square and each 
of E and [^(i) B] has linearly independent rows. 

(c) Suppose F = R with the identity involution. Proceed as in (a), choos- 
ing U, V, and W to be real orthogonal in the two singular value decomposi- 
tions. Thus, A is real orthogonally congruent to a block matrix of the form 
([6]) in which D is square and each of E and [j4(i) B] has linearly independent 
rows. □ 



The regularizing algorithm tells how to construct a sequence of pairs of 
transformations of the square matrices that are sufficient to determine 
the regularizing decomposition of A. Implicit in these transformations is a 
sequence of pairs of *congruences that reduce A in successive stages. After 
the first stage, the *congruences reduce A to the form ([6]). After the second 
stage, if we were to carry out the *congruences we would obtain a matrix of 
the form 





* 





* 


" 




* 


* 


■ 


* 





}m4 











■ 





}m3 


* 


* 


* 


* 


■ 


}m2 














_ 


}mi 



in which the diagonal blocks are square, the * blocks are not necessarily 
zero, and each ■ block has linearly independent rows. Theorem [2] ensures 
that if A is complex, then there are unitary matrices U and V such that 
each of U*AU and AV has the form ([8]), with possibly different values 
for the parameters mi. li A is real, there is a real orthogonal Q such that 
Q'^AQ has the form ®. 

3 * Congruence Invariants and a Reduced Form 

Throughout this section, A E F™^™ and S is a nonsingular matrix. Of 
course, nullity A = nullity S'*j4S', so nullity is a *congruence invariant. The 
relationships 

N{S*AS) = S-^N{A) and N{S*A*S) = S-'^N{A*) (9) 

between the null spaces of A and S*AS, and those of A* and S*A*S, imply 
that 

N{S*A*S) n N{S*AS) = S-^ {N{A*) n N{A)) . (10) 
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We refer to C := dim N{A*) n N{A) as the *normal nullity of A. We let 
:= nullity A, refer lo n := v — C, as the *non-normal nullity of A, and let 
p = m — K — i^. 

It follows from ([9]) and ()10|) that z^, k, and p are *congruence invariants. 
Because 

A 



N{A*)nN{A) = N 



A* 



(11) 



and (" (and hence also k and p) can be computed using elementary row 
operations. 

The parameter mi produced by the regularization algorithm is the nullity 
of A, so it is a *congruence invariant: mi = ly. 

The parameter m2 produced by the regularization algorithm is the rank 
of the block in ([3]). Since has mi columns and full row rank, its 
nullity is mi — m2- Suppose z G F™i and Nz = 0, let y* = [0 z*], and let 
A = SAS* denote the block matrix in ([3]). Then Ay = and y*^ = so 
^ = dim(A^(^*) n A^(^)) = nullity = mi — m2 and hence mi — m2 = C is 
the *normal nullity of A. This means that m-2 = mi — (^ = 1^ — (^ = Kis the 
* non-normal nullity of j4, so m2 is also a * congruence invariant. 

The following lemma ensures that the *congruence class of the square 
matrix ^(x) in Q is also a *congruence invariant. 

Lemma 3. Suppose that a singular square matrix A is *congruent to 



M 





B 





c 


D 


E 











and also to M 



A 



(1) 
C 





B 
D 





E 
0. 



in which D is K-by-K, ID is K-hy-K, and each of E, E, [^(i) B], and [A(^i^ B] 
has linearly independent rows. Then v = v,, n = k, and ^(i) is *congruent 
to A(^i-^, that is, v, k, p, and the *congruence class of the p-by-p matrix ^(i) 
are * congruence invariants of A. 

Proof. The form of M ensures that v is its nullity and that k is its *non- 
normal nullity; is the nullity of M and k is its *non-normal nullity. Since 
M and M are *congruent to A and hence to each other, their nullities and 
*non-normal nullities are the same, so z^ = and k = k. 
Let 



M 



A 



(1) 



C 



B 
D 





E 



and M 



A 



(1) 
C 



B 
D 





E 



US 



[^ij\i,j=i is nonsingular, ^22 is z^-by-z^, and SMS 
SM : 



5*11 


S12 




' M ' 










5*21 


S22 









. S21M 








M, then 



MS- 
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so S21M = 0. Full row rank of M ensures that S21 
and are nonsingular. If we write 511 = [Rij]fj= 
then equating the 1, 2 blocks of SMS* and M_ tells us that 



and hence both 5ii 
^, in which R22 is K-by-K, 



Su 



■ 




R12E 




" ■ 




" 


E 












E 



(5: 



22 J 



which ensures that i?i2 = 0, i?ii and i?22 are nonsingular, and iiiiA(i)i?J^ 
A 



□ 



Lemma El identification of the parameters in the first stage of the regu- 
larization algorithm as *congruence invariants (mi = u and m2 = n), and an 
induction argument ensure that at each stage k = 1,2,... of the algorithm 
the *congruence class of the square matrix j4(fc), the integers m2k-i (the 
nullity of A(fc_i)) and m2k (the *non-normal nullity of ^(fe_i)), the number 
of stages r in the algorithm until it terminates, and the *congruence class 
of the final nonsingular square matrix A(^^^ are all uniquely determined by 
the *congruence class of A. All that remains to be shown is that these 
data determine the regularizing decomposition of A according to the rule in 
Theorem [TJ 

The block matrix ([6]) can be reduced to a more sparse form by *congru- 
ence if m2 > 0: the block E may be taken to be [1^2 0] and the blocks C 
and D may be taken to be zero. To achieve these reductions, is is useful to 
realize that if ^ ^ AS adds linear combinations of a set of columns of A 
with index set a to certain columns, and if the rows of A with index set a 
are all zero, then S*A = A, so 5*^45 = AS. 

Lemma 4. // a singular square matrix A is *congruent to a block matrix 
A of the form in which m2 > 0, D is m2-by-m2, and E has linearly 
independent rows, then it is *congruent to 





B 





Om2 













(12) 



Proof. Since rank E = m2 , there is a nonsingular V such that EV = [Im2 0] • 



For 5 = Im-mi ®V we have 



S*AS = AS 





B 





c 


D 


[Im, 











A'. 
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Then, for 



S 



^m—mi 
^ Irni 



and X 



C D 




J^S = S*A'S has the form m 



□ 



A block matrix of the form (|12|) is said to be a *congruence reduced 
form of A if it is *congruent to A, A(i~^ is square, and [^(i) B] has hnearly 
independent rows. There are four possibihties for the p-hy-p matrix A^i^ in 
a *congruence reduced form of A: 



p = 0: Then A is *congruent to 



A 







"12 





[In., 0] 
Orrti 



Since rank^ = m2 and A"^ = 0, its Jordan Canonical Form contains 
m2 blocks J2 and mi — m2 blocks Ji. But A is similar to its Jordan 
Canonical Form via a permutation similarity, which is a *congruence. 



so J 



[mi-m2] „ j[m2] 



J 2 is the regularizing decomposition for A. 



p> and j4(i) = Op, so = nullity ^(i) = p: A is *congruent to 

Omo B 



A 



-"7713 ^ 

Om2 [Im^O] 
Or 



in which B has full row rank. There is a nonsingular V such that 
BV = [Ims 0], so if we let S = Im^ ®V ® Imi, we have 



Now let S 



S*AS 



S*RS 



0, 








[Im; 0] 
Om2 






[V* 0] 



0. 



and compute 








[^ma 0] 

Om2 







0. 



N. 



Then rank N = 7713 + m2 , rank A^^ 
Canonical Form of is 
regularizing decomposition for A. 



and A^^ = 0, so the Jordan 



which is the 
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p > and ^(1) is nonsingular: Let R denote the block matrix in (|12p. 
let 







'mi 1 



and compute 



S*RS 



^(1) 








X 




[Im, 0] 











in which X = —B*A^*^A(^iy Lemma|4]tells us that S*RS is *congruent 
to ([I2]) with B = 0, that is, to M with 



Jm2 





[Im, 0] 



Since rankM = 1712 and = 0, the regularizing decomposition of A 



is ^(1) e j[ 



"11-^.2] 



J, 



[012] 



/9 > and is singular but nonzero: We address this case in the 
next lemma. 



Lemma 5. Let p > and let A^^) be the p-by-p upper left block in a *con- 
gruence reduced form of A. Let and denote the nullity and 
*non-normal nullity, respectively, o/>l(i), and suppose that > 0. Then 
A is *congruent to 





B' 
















[Lm, 0] 
















Vm, 0] 














0m2 


[Irn, 0] 

















(13) 



in which [A(2) B'] has linearly independent rows. The parameters mi, m2, 
ms, and m4, and the *congruence class of A(2) are *congruence invariants 
of A. 

Proof. Step 1: Lemma H] ensures that there is a nonsingular S such that 



A 



(2) 



B' 

Omi 







[/mi 0] 
Orris 



11 



is a *congruence reduced form of ^(i)- Let /)' denote the size of ^(2)- Let 



5 = 50/, 



m2+mi 



and observe that S*AS has the block form 









" ^(2) 


B' 





Bi 





5M(i)5 


S*B 










07714 


[Irn, 0] 


B2 










[Im, 0] 










0r)i3 


Bs 

















































(14) 



in which 



S*B 



Bi 

B2 

B3 



Step 2: Let M denote the upper left 2-by-3 block of the 5-by-5 block 
matrix in (jl4p . The rows of M are linearly independent, so its columns span 
]pp +m4 ^ ^(jj linear combination of the columns of M to the fourth block 
column of (|14p in order to put zeros in the blocks Bi and B2 ■ Complete this 
column operation to a *congruence by adding the conjugate linear combi- 
nation of rows of M to the fourth block row of (jl4p : this spoils the zeros in 
the first four blocks of the fourth block row. Add linear combinations of the 
fifth block column to the first four block columns in order to re-establish 
the zero blocks there; the fifth block row is zero so completing this column 
operation to a *congruence with a conjugate row operation has no effect. 

We have now achieved a *congruence of A that has the form 



R 



' ^(2) 


B' 
















[In,, 0] 
















B3 
















[Im, 0] 














Orrii 



(15) 



in which B3 has linearly independent rows. 

Step 3: Whenever one has a block matrix like that in ()15p . in which 
some of the superdiagonal blocks below the first block row do not have the 
standard form [/ 0] but nevertheless have linearly independent rows, there 
is a finite sequence of *congruences that restores it to a standard form like 
that in (llSp . For example, -B3 in (jlSp has linearly independent rows, so 
there is a nonsingular V such that B^V = [Img 0]. Right-multiply the 4th 
block column of i? by ^ and left-multiply the 4th block row of the result by 
V*. This restores the standard form of the block in position 3,4 but spoils 
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the [I 0] block in position 4, 5, though it still has linearly independent rows. 
Now right-multiply the fifth block column by a factor that restores it to 
standard form (in this case, the right multiplier is V~* © Imi-m2) then 
left-multiply the fifth block row by the * of that factor. If there are more 
than five block rows, continue this process down the block superdiagonal to 
the block in the last block column, at which point all of the superdiagonal 
blocks below the first block row are restored to standard form since the last 
block row is zero. Of course, this finite sequence of transformations is a 
*congruence of R. □ 



The preceding lemma clarifies the nature of the block S in a *congruence 
reduced form p2|) of A: except for the requirement that [^(i) B] have full 
row rank, B is otherwise arbitrary. 

If there are different involutions on F, the same matrix may have a 
different regularizing decomposition for each involution. For example, take 
F = C and consider 

7'- 

I 1 

If the involution is complex conjugation, then N{A) = N{A*) since A is 
Hermitian, ^ = mi = 1, k = m2 = 0, and p = 1; the regularizing decom- 
position of A is [1] © Ji. However, if the involution is the identity, then 
N{A) n N{A'^) = {0}, ( = nil = 0, K = 1712 = I, and p = 0; the regularizing 
decomposition of A is J2- 



4 The Regularizing Decomposition 

If the block A(^2) in ([13]) is singular, repeat the first two steps in the proof 
of Lemma [5] to reduce A further by *congruence and produce the nullity 
777-5 and *non-normal nullity 7776 of ^(2) and a square matrix ^(s). Then 
perform the process described in Step 3 to restore the standard form of the 
superdiagonal blocks below the first block row. 

Reduction of A to a sparse form that reveals all of its singular structure 
under *congruence can be achieved by repeating the three steps in Lemma 
Oto obtain successively smaller blocks ^{3), ^(4)) •••) ^(r) (with successively 
smaller nullities) in which A(^^^ is the first block that is nonsingular. The 
payoff for our effort in deriving a form more sparse than that produced by 
the *congruences implicit in the regularization algorithm alone, e.g., is 
that it permits us to verify the validity of the regularizing decomposition 
asserted in Theorem [TJ 
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Theorem 6 (Regularizing Decomposition). Let A be a given square sin- 
gular matrix over F. Perform the regularization algorithm on A and ob- 
tain the integers T,mi,m2, ■■■,m2T and a nonsingular matrix ^(r)- Then 
r, mi, 777,2, and the *congruence class of A^^-j are *congruence invari- 
ants of A. Moreover, 

(a) (Canonical sparse form) A is *congruent to A^^^ © N , in which 



N 



0, 



m2T 



\Im2T 

0. 



0] 



m2T-i 



0, 



0] 



m2 



[Im2 0] 



(16) 



has all of its nonzero blocks [/m2r 0]) • • • ) [-^m2 0] first block superdiag- 

onal, and each block [Imk 0] Tife-fey-mfe-i; A; = 2, 3, ... , 2r. 
(b) (Existence) A is *congruent to A^^^ © M , in which 



jYiri'ZT-l 



-m2T 



J. 



[m2r] 



2t 



(17) 



(^cj (Uniqueness) Suppose A is *congruent to B®C , in which B is nonsingu- 
lar and C is a direct sum of nilpotent Jordan blocks. Then B is *congruent 
to yl(^) and some permutation of the direct summands of C gives M . 
(d) (Unitarily reduced form) Consider the three cases (a) F = C and the 
involution is complex conjugation, or ((3) F = C and the involution is the 
identity, or (j) F = M and the involution is the identity. Then depending 
on the case there is (a) a complex unitary U , or ((3) a complex unitary V , 
or (j) a real orthogonal Q such that (a) U*AU or (f3) V'^AV or (7) Q'^AQ 
has the form. 



}m2T 
}m2T-i 



B2T+1 








* 


* 


B2t 












* 




Bj 







Be 




* 




B5 



* B4 




* 




B. 



* B2 




(18) 



}m4 

}"^3 
}m2 
}mi 
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in which all 2r + 1 diagonal blocks B2T+1, *, 0, . . . , *, 0, *, are square, B2T+1 
is nonsingular, and each of B2, ■ ■ ■ , has linearly independent rows. The 
integers T,mi, . . . ,m2T o,re the same as those in U6\) and ( [j7| j. The equiv- 
alence class (under complex *congruence, complex T -congruence, or real 
T-congruence, respectively) of B2T+1 is the same as that of A^^^y In the 
principal suhmatrix of hl^) obtained by deleting the block row and column 
containing B2T+1, replacing all blocks denoted by stars with zero blocks and 
replacing each Bi with [I 0] produces the matrix N in [Tt 



Proof. The *congruence invariance of the parameters and r, as well as 
the *congruence class of A(t-) have already been established. The form of 
N is the outcome of repeating the reduction described in Lemma [5] until 
it terminates with a block A(^^-^ that is nonsingular. The only issue is the 
explicit description of the Jordan block structure in (|17|) . 
Notice that 

[Im, 0] 0] = [Im, 0] 



"ife_i 



and hence 



mfc-2 



"lfe-2 



[/m^.-i 0] 



V3 [Ira, 0] 



OjTll 



has its nonzero blocks [Im2T ^]-, ■ ■ ■ A^ma 0] in the second block super- 
diagonal. In general, is a 0-1 matrix that has its nonzero blocks 
[Im2T 0]) • • • ) [-^m,fc+i 0] in the kth block superdiagonal. The structure of 
the powers A^'^ ensures that the rank of each is equal to the number of its 
nonzero entries, so 



rank N ' = m^+i + • • • + ^27- 



1, 



,2r- 1 



(19) 



and N'^'^ = 0. The list of multiplicities of the nilpotent Jordan blocks in 
the Jordan Canonical Form of (arranged in order of increasing size) is 
given by the sequence of second differences of the sequence {rankiV'^}^^^ 
[5l Exercise, p. 127], which is mi — m2, m2 — m^, — m^, etc. The direct 
sum of nilpotent Jordan blocks in (I17p is therefore the Jordan Canonical 
Form of N. 
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The final step in proving (jl7p is to show that the Jordan Canonical Form 
of N can be achieved via a permutation similarity, which is a *congruence. 
A conceptual way to do this is to show that the directed graphs of the two 
matrices M and N are isomorphic. 

The directed graph of is a linear chain with k nodes Pi,. . . ,Pk in 
which there is an arc from Pi to P^.+i for each i = 1, . . . ,k — 1, so the 
directed graph of M is a disjoint union of such linear chains. There are 
iTT'k ~ i^k+i chains with k nodes for each k = 1, . . . , 2t. 

To understand the directed graph of one can begin with any node 
corresponding to any row in the first block row. Each of these m2r nodes 
is the first in a linear chain with 2t nodes. In the second block row, the 
nodes corresponding to the first rows are members of the linear chains 
associated with the first block row, but the nodes corresponding to the last 
m2r-i — "i2r TOWS begin their own linear chains, each with 2t — 1 nodes. 
Proceeding in this way downward through the block rows of N we identify 
a set of disjoint linear chains that is identical to the set of disjoint linear 
chains associated with M. A permutation of labels of nodes that identifies 
the directed graphs of M and N gives a permutation matrix that achieves 
the desired permutation similarity between M and N. 

The uniqueness assertion follows from (a) our identification of all the 
relevant parameters as *congruence invariants of A and (b) uniqueness of 
the Jordan Canonical Form. 

Finally, the assertions about the unitarily reduced form (jlSp follow from 
the regularizing algorithm in Section [2] and the proof of Theorem [2j When 
the regularizing algorithm is carried out with unitary transformations, the 
result is a matrix of the form (1181) . of which ([8]) is a special case. □ 



5 Regular izat ion of a *Selfadjoint Pencil 

Theorem [6] implies that every *selfadjoint matrix pencil A + XA* has a reg- 
ularizing decomposition ([2]) with a *selfadjoint regular part. The algorithm 
in Section [2] can be used to construct the regularizing decomposition, and 
if F = C with either the identity or complex conjugation as the involution 
(respectively, F = M with the identity involution), the construction can be 
carried out using only unitary (respectively, real orthogonal) transforma- 
tions. We emphasize that the involution on F may be the identity, so the 
assertions in the following theorem are valid for matrix pencils of the form 
A + XA^. 
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Theorem 7. Let A + XA* be a *selfadjoint matrix pencil over F and let A 
be *congruent to A(^^^ ©M, in which ^(^-j is nonsingular and M is the direct 
sum of nilpotent Jordan blocks in |j7| j. Then there is a nonsingular S such 

K and 



® {J2r + XJlf"^'^ 



(20) 



that S{A + XA*)S* = {A(^r) + ^^1^)) 

Moreover, each singular block + A jj' may be replaced by 

{ [Fi + XGi) © {Gj + XFf) if k = 21 -I is odd 
y [Ji + Xli) © {h + XJi) ifk = 2e is even. 

Use of the blocks ()20p instead of the corresponding Jordan blocks is 
justified by the following lemma. 

Lemma 8. Jfc + A jj is strictly equivalent to i20\). 

Proof. If there is a permutation matrix S such that 



SJkS 



Mp 





Ff, 



gJ 



if /c = 2^ - 1 is odd 



\i k = 21 is even, 





h 

Jrn 

then S{Jk + AJj)5^ is strictly equivalent to (pO]) . To prove the existence 
of such an S, we need to prove that Mi and A^^ can be obtained from Jk 
by simultaneous permutations of rows and columns, that is, there exists a 
permutation / on {1, 2, . . . ,k} that transforms the positions 

(1,2), (2,3), ...,{k-l,k) 

of the unit entries in to the positions 

(/(l),/(2)), (/(2),/(3)), {f{k-l),f(k)) (21) 

of the unit entries in Mi if k = 21 — 1 or in Ni ii k = 2i. To obtain the 
sequence ([2T]) , we arrange the indices of the units in 



Mp 








" 





1 •. 










1 


1 







10 







{{2i-l)-hy-{2£-l)) 
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as follows: 

{£,2i-l), {2£-l,£-l), {e-l,2£-2), {2i-2,i-2), 

(2,1 + 1), (£+1,1), 

and the indices of the units in iVj. as follows: 

(1,1 + 1), (i + 1,2), (2,e + 2), (£ + 2,3),...,(2£-l,e), (£,2i). 

□ 
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